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Abstract 

We review the currently known universality classes of continuous phase transitions 
to absorbing states in nonequilibrium systems and present results of simulations 
and arguments to show how the blockades introduced by different particle species 
in one dimension cause new robust classes. Results of investigations on the dynamic 
scaling behavior of some bosonic spreading models are reported. 
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1 Introduction 

Critical universality is an attractive feature in statistical physics because wide 
range of models can be classified purely in terms of their collective behav- 
ior into classes. While in equilibrium systems the factors determining these 
classes are quite well understood, in nonequilibrium systems the situation is 
less clear. The 1+1 dimensional reaction-diffusion systems are important for 
the understanding of the universality classes of non-equilibrium system be- 
cause order-disorder phase transitions may occur in such a low dimension 
unlike in equilibrium systems [1]. However, for the ordered state to be sta- 
ble it must exhibit only small fluctuations. All of this kind of transitions are 
such that the ordered state is "absorbing", i.e. when the system falls into it, 
it can not escape. In fact for a long time only such phase transitions have 
been known where the absorbing state is completely frozen. A few universal- 
ity classes of this kind are known [2] , [3] , the most prominent and the first one 
that was discovered is that of the directed percolation (DP) [4]. An early hy- 
pothesis [5] was confirmed by all examples up to now. This claims that in one 
component systems exhibiting continuous phase transitions to single absorb- 
ing states (without extra symmetry, inhomogeneity or disorder) short ranged 
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Table 1 



Summary of known 1+1 D universality classes. DCF, NDCF denote models with 
coupled (non) diffusive conserved fields. 



CLASS ID 


features 


degree of knowledge 


references 


DP 


time reversal symmetry 


RG to e^, series exp. 


[1-61 


PCP 


broken time reversal symm. 


RG, simulations 


[81 


NDCF 


trlobal conservation 


RG simulations 


[10-12] 


CDP 


Comnart DP Olanhpr Tsinp r 


pxactlv soIvpH 


[13-15] 


PC 


naritv pnrm -U Tjn wrmnptrv 


T?(~J si mnlflt.invm 

j. i vi . oiiiiuicx uiuiio 


[16-24] 


N-BARW2 


N-comp. parity conservation 


RG, simu., DMRG 


[18,25,26] 


PCPD 


binary production 


simulations 


[27-33] 


DCF 


global conservation 


RG, simulations 


[34-37] 


N-BARW2s 


N-comp. sym. BARW + exc. 


simulations, MF 


[25,38-40] 


N-BARW2a 


N-comp. asym. BARW + exc. 


simulations, MF 


[38,25] 



interactions can generate DP class transition only. The Lagragian of the field 
theory of the DP process L{<p{x, t),ip{x, t)) exhibits a time reversal symmetry 
(j>(x, t) <-> ip(x, —t). This results in relations among scaling exponents (see [3]). 



The same static scaling behavior was observed in frozen multi-absorbing state 
systems as well, like in the pair contact process (PCP) [8]. Here the time re- 
versal symmetry is broken owing to an extra term in the Lagrangian describ- 
ing the effect of frozen particles generating long time memory in the system 
[9]. Besides DP a few more universality classes have been established in the 
last decade. In models exhibiting particle annihilation fluctuating absorbing 
states with a single wandering particle can occur [7]. In Table 1 we have tried 
to summarize the most well known absorbing state phase transition classes 
of disorder- free, homogeneous models with short range interactions. Those, 
which are below the horizontal line exhibit fluctuating absorbing states. For 
more details of field theoretical symmetries and their relation to hyperscaling 
laws see [9]. 

A major problem of these models is that they are usually far from the critical 
dimension and critical fluctuations prohibit mean-field (MF) like behavior. 
Further complication is that bosonic field theoretical methods can not de- 
scribe particle exclusion that may obviously happen in ID. The success of the 
application of bosonic field theory for models shown in the first part of the 
table is the consequence of the asymptotically low density of particles near the 
critical point. However in multi-component systems, where the exchange be- 
tween different types is non-trivial, bosonic field theoretical descriptions may 
fail. Also in case of the binary production (PCPD) models [27] it predicts a 
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non-vanishing density at the transition point and diverges in the active phase 
contrary to the lattice model version of hard-core particles [28-33] . Fermionic 
field theories on the other hand have the disadvantage that they are non- 
local and results exist for very simple reaction-diffusion systems only [41,42]. 
Other techniques like independent interval approximation [43], empty inter- 
val method [44] or density matrix renormalization (DMRG) [45] are currently 
being developed to be able to solve mult i- component reaction diffusion mod- 
els. The main aim of this work is to give an overview and present simulation 
results for these systems. 



2 Critical dynamical behavior of single-species spreading processes 

The bosonic field theory of the directed percolation was established by [46]. 
Perturbative e = 4 — d renormaliztion group analysis [5] up to two-loop order 
resulted in estimates for the critical exponents. For the density decay exponent 
it gives 

a B = 1 - e /4 - 0.01283e 2 (1) 

which results in a B = 0.13453 in one dimension. This differs from the most 
precise simulation and series expansion results of DP a = 0.1595(1) [49] with 
about 18%. An attempt for fermionic field theoretical solution in Id (that does 
not allow multiple occupancy of sites) was shown in [42] but run into severe 
convergence problems and has not resulted in precise quantitative estimates 
for critical exponents. Although the bosonic field theory is expected to be valid 
owing to the asymptotically low density at criticality it has never been proven 
rigorously. Contrary, a recent study on a special one-dimensional DP model 
predicts different bosonic and fermionic field theoretical results [47]. As far as 
we know all simulations of directed percolation (and other spreading models) 
have been done on lattices allowing single occupancy. Moreover there are no 
more precise estimates for the critical exponents of the Id bosonic model than 
those of the e 2 expansion cited above. Therefore we decided to perform some 
simulations with bosonic particles to check the scaling behavior at criticality. 

First we investigated the following one-dimensional branching and annihilating 
process with one offspring (BARW1) 

A^2A 2A^(D A%^ X %A (2) 

in such a way that the branching and the coagulation processes happen in 
place. The simulations were run on L = 10 6 systems with periodic boundary 
conditions and with initially randomly distributed 0-s and As of probability 
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0.142 




Fig. 1. Local slopes of the density decay in a bosonic DP model. Different curves 
correspond to A = 0.12883, 0.12882, 0.12881 0.1288, 0.12879 (from bottom to top). 

1/2. The density of particles (p(t)) is measured up to t max = 2 x 10 5 Monte 
Carlo steps (MCS) and averaged over 10 4 samples. By fixing the branching 
rate to o = 0.1 and varying the annihilating rate we determined the critical 
point A c with the local slope analysis of data 

-\n[p(t)/p(t/m)] 
Qeff{t) ~ Umj (3) 



(where we use m = 8 usually). In the t — > oo limit the critical curve goes 
to exponent a by a straight line, while in sub(super)-critical cases they veer 
down(up) respectively. As Figure 1 shows the critical point is at A c = 0.12882(1) 
with a = 0.165(5) that agrees well with other simulation and series expansion 
results 0.1595(1) [49]. By simulating on smaller lattice sizes (L = 20000) we 
found strong finite size effects. 

Next we performed the same analysis for the Id two-offspring version BARW 
(BARW2) 

A^3A 2A^(& A0 1_ A _A 0A (4) 



that exhibits mod 2 particle number conservation and its critical behavior is 
known to belong to the PC class [16,18]. Now we fixed A = 0.2 and as Fig. 2 
shows the critical point is at a c = 0.04685(5) with the corresponding decay 
exponent a = 0.290(3). This value agrees with that of the PC class again 
0.285(2) [17]. 
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Fig. 2. Local slopes of the density decay in a bosonic BARW2 model. Different 
curves correspond to a = 0.466, 0.468 0.469, 0.47 (from bottom to top). 

Finally we performed bosonic simulations for the sub-critical behavior of the 
following Id binary spreading process 

2A A 3A 2A^(& A(& 1_ ^ A (&A (5) 

In the fermionic version of this model continuous phase transition was found 
[28-30] with new critical behavior. On the other hand the field theory of 
the bosonic model [27] predicted discontinuous phase transition at 2A = a 
with diverging particle density in the active phase. Our present simulations 
confirm this. While in the inactive phase the 2 A — > process (see Section 3.1) 
governs the evolution it is guessed [48] that at the transition point the 3A — > 
dominates with the following scaling law in Id 

p(t)oc(ln(t)/t) 1 / 2 + 0(lA) . (6) 

The simulations at o = 0.2 in systems with size L = 10 6 confirmed these 
expectations as shown on Fig. 3. At the transition point for t > 4 x 10 4 MCS 
we applied the form eq.(6) and found a good fitting with pit) = ((0.0266(1) — 
0.1765(1) ln(t)) /t) 1 / 2 . Though the identification of a universality class requires 
the determination of three independent critical exponents we believe that the 
above numerical data give enough support for the expected critical behaviors. 
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Fig. 3. Density decay in the bosonic annihilation-fission model. The upper curve 
corresponds to the inactive phase A = 0.5, the lower one to the transition point 
A = 0.1. The dashed line shows logarithmic fitting of the form eq.(6). 

3 Critical dynamical behavior of multi-species annihilating random 
walks 

First we recall some well known results (see refs. in [50]) for simple reaction- 
diffusion systems, that lack particle creation. From the viewpoint of phase 
transitions this describes the behavior in the inactive phase or as we shall 
show in the N-BARW models right at the critical point. 

3.1 Models A + and A + B ^0 

The simplest reaction-diffusion model A + A — > (ARW) - in which identical 
particles follow random walk and annihilate on contact is exactly solvable in 
ID [51,52], the particle density decays as 

p(t)<xrV* (7) 

ARW was also shown to be equivalent to the A + A — > A model by [53] and 
renormalization group approach provided universal decay amplitudes to all 
orders in epsilon expansion. It was also shown [14] that the motion of kinks in 
the compact version of directed percolation (CDP) [13] and in the Glauber- 
Ising model [15] at zero temperature are also described exactly by (7). 
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Two types of particles undergo diffusive random walk and react upon contact 
to form an inert particle in the simplest two-component reaction-diffusion 
model A + B — > 0. For d < 4 and for equal initial density of A and B particles 
(po), which are randomly situated in space the density decays asymptotically 
as [54] pA(t) = PB(t) oc Cd^/p^t" 1 ^ where Cd is a dimension dependent con- 
stant. This slow decay is due to the formation of clusters of like particles that 
do not react and their asymptotical segregation for d < 4. The asymptotically 
dominant process is the diffusive decay of the fluctuations in the initial con- 
ditions. Since this is a short ranged process the system will have a long-time 
memory - appearing in the amplitude dependence - for the initial density. 



3.2 The effect of exclusion 

A simple reaction-diffusion model of two types A + A ^ B + B ^ $ with 
exclusion AB <A BA is offered by the Generalized Domany-Kinzel (GDK) 
stochastic cellular automaton model of Hinrichsen [24] introduced originally 
to realize phase transitions from active (Ac) to multiple inactive absorbing 
states: II, 12. We investigated numerically [55] this model in a special point 
of its phase diagram where compact domains of II and 12 grow separated by 
Acll = A and AcI2 = B kinks that can not penetrate each other (CDP2). 
There are special pairwise initial conditions in the model (because the domains 
are bounded by kinks of the same type): ....A...A...B.B..B B..A..A.. 

The dynamical behavior of critical systems is usually investigated from two ex- 
treme initial conditions: (a) homogeneous system with randomly distributed 
species (b) empty system with a single initial seed of particles. In the for- 
mer case (a) and pairwise initial conditions our numerical simulations show 
a density decay of kinks (or particles of the corresponding reaction-diffusion 
system) p oc t~ a characterized by a power-law with an exponent larger than 
a = 0.5 that would have been expected in case of two copies of ARW sys- 
tems that do not exclude each other. Furthermore the deviation of a from 
0.5 showed an initial density dependence. We provided a possible explanation 
based on symmetry between types that a marginal perturbation emerges here 
that causes this non-universal scaling. We showed an analogy with the DP con- 
fined by parabolic boundary conditions [56] by assuming that the kinks exert 
a parabolic space-time confinement on the decaying domains. Non-universal 
scaling can also be observed at surface critical phenomena similarly to here 
where kinks produce 'multi surfaces' in the bulk. However simulations and 
independent interval approximations on a similar model predict a t' 1 ^ 2 / ln(t) 
behavior [57]. Perhaps a fermionic field theoretical study could help to under- 
stand better the situation. In case of random distribution of As and Bs the 
density decays slower owing to the accumulation of AB pairs. An exact map- 
ping onto the A + B — > model predicts Pa = Pb oc t^ 1 ^ that was confirmed 
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by simulations [25]. 



In case (b) above one usually measures the survival probability of clusters 
originated from a seed. In critical systems this scales like P(t) oc t~ 5 . When 
we inserted a seed of /2-s in the see of II and Acs of the GDK model at the 
CDP2 point we found very strong dependence of 5 on po(Il) [55]. Also the 
above picture with parabolic boundary conditions has gained further support 
here. 



4 Phase transition generated by offspring production 

If we add particle creation by branching to the systems mentioned before we 
may expect a phase transition, where a steady state with constant particle 
density occurs. For iV component parity conserving systems (N-BARW2, N = 
2 types, two offsprings) it was shown by field theory [18] that the A — > 3A 
like processes are not relevant and the models with A — > A2B and B — > B2A 
(ordering of offsprings is not relevant) branching terms do exhibit continuous 
phase transitions at a = branching rate. The universality class is expected to 
be independent from N and to coincide with that of the N — > oo (N-BARW2) 
model. The critical dimension is d c = 2 and for d — 1 the exponents are 

(3 = 1, Z = 2, a = 1/2, u\\=2, u ± = l (8) 

exactly. 

4-1 The parity conserving 2-BARW2 model with exclusion 

The effect of exclusion in the 2-BARW2 model was investigated in [25] and 
for d = 2 the field theoretical predictions were confirmed. In one dimension 
however two phase transitions of different types were observed depending on 
the arrangement of offsprings relative to the parent. Namely if the parent 
separates the offsprings: A —> BAB (2-BARW2s) the steady state density will 
be higher than in the case when they are created on the same site: A — > ABB 
(2-BARW2a) for a given branching rate because in the former case they are 
unable to annihilate with each other. This results in different order parameter 
exponents for the symmetric and the asymmetric cases (f3 s = 1/2 and j3 a = 2). 
This is in contrast to the widespread beliefs that bosonic field theory (where 
AB <-> BA is possible) can well describe these systems because in that case 
[18] the scaling behavior is different (8). This finding led [38] to the conjecture 
that in one-dimensional reaction-diffusion systems a series of new universality 
classes should appear if particle exclusion is present. Note however that only 
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7 


OL 


R 


N-BARW2 


2 


2 


1/2 


1 


N-BARW2s 


2.0(1)|0.915(2) 


4.0(2)|1.82(2)* 


0.25(1)|0.55(1)* 


0.50(1) 


N-BARW2a 


8.0(4)|3.66(2) 


4.0(2)|1.82(2)* 


0.25(1)|0.55(1)* 


2.05(10) 



Table 2 

Summary of critical exponents in one dimension in N-BARW like models. The non- 
blocking data are quoted from [18]. Data divided by "|" correspond to random vs. 
pairwise initial condition cases [55]. Exponents denoted by * exhibit slight initial 
density dependence. 

the off-critical exponents are different, since in both cases the transition is at 
o" = and the on-critical ones are those described in Sect. (3.2). In [25] a set 
of critical exponents satisfying scaling relations have been determined for this 
two new classes shown in Table 2. 



4-2 The parity non- conserving model with exclusion (2-BARW1) 



Hard-core interactions in the two-component, one-offspring production model 
(2-BARW1) were investigated in [39]. Without interaction between different 
species one would expect DP class transition. By introducing an AB <A BA 
blocking in the 



A AA B ^ BB (9) 

AA ^ BB ^ (10) 

model a DP class transition at a = 0.81107 was been found indeed. On the 
other hand if we couple the two sub-systems by the production 



A^AB A^BA (11) 
B ^ AB B^BA (12) 
AA ^> BB (13) 



without exclusion bosonic field theory [58] predicts unidirectionally coupled 
DP class [59] transitions. With hard-core exclusion a continuous phase tran- 
sition will emerge at rate a = - therefore the on-critical exponents will be 
the same as described in Sect. (3.2) - and the order parameter exponent was 
found to be (3 = 1/2. This assures that this transition belongs to the same 
class as that of the 2-BARW2s model. The parity conservation law that is 
relevant in case of one component systems (PC versus DP class) is irrelevant 
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here. This finding reduces the expectations for a whole new series of universal- 
ity classes in ID systems with exclusions. In fact the blockades introduced by 
exclusions generate robust classes. In [39] a hypothesis was set up that in cou- 
pled branching and annihilating random walk systems of N-types of excluding 
particles for continuous transitions at a = two universality classes exist, 
those of 2-BARW2s and 2-BARW2a, depending on whether the reactants can 
immediately annihilate (i.e. when similar particles are not separated by other 
type(s) of particle(s)) or not. Recent investigations in similar models [40,60] 
are in agreement with this hypothesis and the extension to binary spreading 
processes was proposed in [61]. 



5 Summary 



We have shown how blockades generated by particle exclusion in some one- 
dimensional reaction-diffusion type systems affect critical scaling behavior. 

We confirmed by bosonic simulations that the single component unary BARW1 
and BARW2 processes exhibit the same dynamical scaling as the correspond- 
ing models with exclusion. In case of the one component binary spreading 
process the bosonic and fermionic versions differ. While the former one has a 
discontinuous phase transition the latter exhibits a continuous transition. We 
determined numerically the dynamic scaling behavior of the bosonic version. 

Without particle production initial condition dependent dynamical scaling was 
found in multicomponent ARW models. This critical behavior was shown to 
give the set of on-critical exponents if we generate phase transition by particle 
branching since the transitions occur at the a = rate. The off-critical expo- 
nents were found to be insensitive to parity conservation law and depend only 
on the spatial arrangement of parent and offspring. Two generic universality 
classes: 2-BARW2a and 2-BARW2s were explored by numerical simulations. 
In these models the hard-core exclusions generate extra fluctuations in the 
absorbing state and the passive steady states posses long-range correlations. 
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